Radiative damping of surface plasmon resonance in spheroidal metallic nanoparticle 

embedded in a dielectric medium 
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The local field approach and kinetic equation method is applied to calculate the surface plasmon 
radiative damping in a spheroidal metal nanoparticle embedded in any dielectric media. The radia- 
tive damping of the surface plasmon resonance as a function of the particle radius, shape, dielectric 
constant of the surrounding medium and the light frequency is studied in detail. It is found that 
the radiative damping grows quadratically with the particle radius and oscillates with altering both 
the particle size and the dielectric constant of a surrounding medium. Much attention is paid to 
the electron surface-scattering contribution to the plasmon decay. All calculations of the radiative 
damping are illustrated by examples on the Au and Na nanoparticles. 
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I. INTRODUCTION 

The surface plasmons (SPs) excited by an electro- 
magnetic radiation falling on a metal nanoparticle (MN) 
are still of great fundamental interest ^^"^ since their pro- 
nounced local resonances, whose position, shape and in- 
tensity can be tuned over wide spectral range by varying 
the size and shape of the MNs or by changing the sur- 
rounding medium. 

Once excited, plasmon oscillations can damp non- 
radiatively by absorption caused by electron-phonon in- 
teractions, and/or radiatively by the resonant scatter- 
ing process. ^'^ The electron-phonon interaction with de- 
creasing MN size becomes more and more ineffective 
due to kinematical restrictions imposed on the energy- 
momentum conservation laws.^ Nevertheless, the line 
broadening has been observed in experiments with MNs 
scattering and absorptiouj^i^r— This means that besides 
electron-phonon interaction the other damping mecha- 
nisms should be studied as well in order to interpret the 
observed SP line broadening. 

Usually, both the surface and the radiative damping 
mechanisms play an important role in the SP decay. In 
Ref. [13 it was shown that both the electron surface scat- 
tering and radiative damping can make significant con- 
tribution to the linewidth T (the full width at a half of 
maximum of the surface plasmon resonance (SPR)). 

In the MNs of a smaller radii, the penetration depth 
of the plasmon field reduces and becomes more localized 
near the surface.^ This is due to the fact that the electron 
screening is increased as the particle radius is reduced. 
As a result, the bulk-induced loss processes play only a 
minor role and the electronic excitations generated by the 
surface potential dominate. On the other hand, if elec- 
trons are confined in nanoparticles, where the mean free 
path (MFP) of the electrons becomes comparable with 
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the size of the particles, the electron-surface scattering 
comes into play^ and the surface acts as an additional 
scatterer. 

The radiative damping of the SPR is an important 
parameter since it would help to analyze the specificity 
of the transformation of the collective electron oscillation 
energy into the optical far field. 

The interest to treat of the SPR and radiative damp- 
ing is maintained as well by the investigation of the local 
field enhancement, an effect which increases the inten- 
sity of the incident light near the MNs surface by sev- 
eral orders of magnitude.— 1^ Since a lot of devices in- 
corporating MNs gains from this effect, T is treated as 
a main parameter in applications such as field concen- 
tration for nanopatterning with nanowireS)^ plasmonic 
nanolithography,— and near-field optical microscopy,— 
astigmatic optical tweezers,— surface enhanced Raman 
scatteringi^ etc. 

The damping of SPR in MNs previously has been well 
studied theoretically.^-"^^ The influence of a nearby sur- 
face on the plasmon resonance of a metallic nanoparticle 
of finite size has been studied in the work [^^ with an 
account for the effects of both radiative and evanescent 
surface-reflected waves. In Ref. Hi] the oscillations of 
the plasmon linewidth as a function of the radius of the 
nanoparticles were obtained from numerical calculations 
based on the time dependent local density approximation 
without regard for radiative damping. The temperature 
effect on the radiative lifetime of the surface plasmon 
was studied in the work [22l. To evaluate size and shape 
dependent dielectric functions and extinction spectra for 
MNs with various shape, the different empirical formulas 
were proposed for electron MFPi^ 

The radiative damping calculated in Refs. [2^ - 1261 is 
proportional to the product of the polarizability (propor- 
tional to particle volume) times = {2tt/X)^, where k is 
wave number and A is the wavelength. The calculations 
have been conducted in the electrostatic approximation, 
assuming V/X'^ ^ 1, where V is the MN volume. 

In the works mentioned above was dropped from ac- 
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count the fact that the role of an external electric field 
(which causes the dipole oscillations of electrons and the 
radiative damping) can play the inner electric field in 
MN. Besides, in these works the surface and radiative 
decay have not been considered simultaneously for non- 
spherical MN. 

The purpose of the present paper is to calculate the 
radiative damping of the surface plasmon under the ef- 
fect of inner electric field, for the case when the MFP 
of the electrons is larger than the particle size and their 
scattering on the particle surface plays an important role. 

Within the framework of a local field approach and 
the kinetic equation method, we have found that the 
linewidth of SPR increases quadratically with the parti- 
cle radius and exhibits oscillations as a function of the 
nanoparticle radius and the refractive index of a sur- 
rounding medium. 

The rest of the paper is organized as follows. In Section 
II the theoretical background to the problem is presented. 
Section III is devoted to the study of radiative damping. 
Section IV contains the calculation of the conductivity 
tensor and the discussion of obtained results. Section V 
contains the summary and conclusions. 



II. SURFACE PLASMON LINEWIDTH 

An interaction of light with a MN embedded in a 
medium is studied in the framework of classical optics, 
assuming that the particle and the medium are contin- 
uous, homogeneous, and characterized by their dielec- 
tric function. To overcome the problem connected with 
an inhomogeneous line broadening (due to the size and 
shape distribution within the particle ensemble), we re- 
strict ourselves only to a single MN, which directly yields 
the homogeneous linewidth. ^'-^^ 

In general, the linewidth can be decomposed into con- 
tributions from the bulk dielectric constant, surface scat- 
tering, and radiative damping:"' F = Fb -I- Fg -I- Fiad- 

In the classical case of free electrons in bulk metal, 
the damping Fb(= v) is due to the inelastic scattering of 
the electrons on phonons, lattice defects, or impurities {v 
refers to the electron collision frequency), which shorten 
the MFP.^^ In this case, the relation Fb — vp/loo holds, 
where vp the Fermi velocity and loo is the MFP of 
conduction electrons in the bulk (when the MN radius 
tends to infinity). 

Similarly, to estimate the surface effect in electron- 
surface scattering, the following empirical relatior^^^ 



Ictf 



(1) 



is often used, where A is a phenomenological factor and 
IcS is a reduced effective MFP. For the sphere, the clas- 
sical theory gives IcB = R for an isotropic scattering or 
IcS — 4_R/3 for a diffusive scattering, while the quantum 
theory in a box model yields Ics — 1.16i? or l^s = 1.33i?, 
where R is the radius of a spherical MNi^ 



The bulk scattering does not cause appreciable attenu- 
ation for the high frequency case, when wr 3> 1, where r 
is the relaxation time. On the other hand, the contribu- 
tion to damping from surface scattering is negligible only 
if Zoo/2i? ^ Vl + w^T^. Thus, the surface contribution 
is important if ujt 3> max(l, lao/'2R). 

The radiative contribution can be estimated with the 
help of the following expression^ 



F,ad = 2hkV, 



(2) 



where V is the nanoparticle volume and k is an another 
phenomenological constant to be taken from the experi- 
mental data. 

To account both the surface and the radiative ef- 
fects for a MN embedded in a homogeneous, transparent 
medium, the dielectric permittivitj^i^ 



e'{uj) = ei„tcr(w) + (l - 



(TVfTl+r;:!? 



(3) 



is usually used, where e' and e" are, respectively, the real 
and imaginary parts of the dielectric function of the par- 
ticle material, eintcr accounts for the mterband electron 
transitions, = Ai:nee^ /m, and Ue refers to the elec- 
tron concentration. The expression in the parentheses is 
applied for the mfraband electron transitions. ^-S 

One can see that the effect of the surface is reduced 
simply to the addition of a term 7s = ^s{lcs) in the 
denominator of Eq. ([3]) in the form of Eq. ([T]), and the 
radiative effect is accounted simply by the term 7^ in the 
same denominator. Formulae ([ij and ([2]) can be applied 
only to the MNs of a spherical shape in the case when 
the MFP of electrons I is smaller than the particle size 
d. If the shape of MN differs from the spherical one or 
the inequality I > 2R takes place, then the expressions 
(P) and ^ can no longer be used. 

We will consider the MNs with a moderate sizes for 
whose, on the one hand, the condition I > 2R is still 
fulfilled (the collisions of the conduction electrons with 
the particle surface remains to be an important relax- 
ation process), and on the other hand, with the sizes 
enough large to account for the dissipation of the elec- 
tron energy due to the emission of electromagnetic waves 
by plasmons, so called radiation damping. 



III. RADIATIVE DAMPING 

The problem of a damping of the electron energy due 
to the radiation of a portion of the collective electron 
oscillation energy into the optical far field has been ex- 
tensively studied in the literaturej^iiii^^i^^ 

As we can see from Eq. ([T]) the surface damping de- 
pends on a particle size. The relative contributions from 
the radiative damping through the resonant scattering 
and absorption also strongly depend on the particle size. 
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In particular, it is known^i^ that the plasmon absorption 
is the only process taken place in small particles, whereas 
both the absorption and the scattering are present in 
large particles, with the latter becoming more dominant 
as the particle size increases. The phenomenon is based 
on an interplay of the dissipative and radiative damping. 

An another assumption that the wavelength of the ab- 
sorbing light A is far above the characteristic size of the 
nanoparticle (about 25 nm for gold particles^) allows us 
to treat the MN as being immersed in a spatial uniform, 
but oscillating in time electric field. This implies that 
EM field around the MN can be considered as homoge- 
neous and across the particle as uniform, such that all 
the conduction electrons move in-phase producing only 
dipole-type oscillations. 

The external electric field E*^"-* exp {—iuit) induces an 
inner (potential) electric field Ein inside the particle 
which is coordinate independent. The field Ejn can be 
linearly expressed in terms of E^^' by employing the de- 
polarization tensor. In terms of principal axes, the depo- 
larization tensor, which coincide with with the principal 
axes of the ellipsoid, the relation between components of 
external and inner fields look as^ 



E 



(0) 



(4) 



Ej^in are the components of the electric field inside the 
MN and Lj are the principal value of the j-th compo- 
nent of the depolarization tensor that is also known as 
a geometric factors. The explicit expressions of Lj for 
a MN with a particular shape can be found elesewhere 
(see, e.g., Refs. [soj, [ij, and [Slj). The complex dielec- 
tric permittivity of the particle material is denoted by 
e(= e' + e") and £,„ refers to the dielectric constant of 
the surrounding medium. 

Electrons are accelerated in the presence of an electric 
field inside the MN. It is well known from the classical 
electrodynamics^ that accelerated charges emit electro- 
magnetic radiation in all directions. To calculate the line 
broadening that is entirely caused by an increase of Trad 
due to the radiative effect, we will use the time depen- 
dence of a classical dipole oscillator. The force of a decel- 
erative radiation of a dipole under an inner electric field 
(Eq. ([3])) can be represented as 



■ rad 



it) 
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3?' 



[l + L(e/e„,-l)] d(t), (5) 



where d is the MN dipole moment. The minus sign 
means that this force is opposit to the dipole moment 
direction. In the case of one electron and a medium 
with e,„ = £ = 1, Eq. ([5]) transforms into the well- 
known expressions from the classical electrodynamics;^ 
The linewidth due to the radiative damping of dipole vi- 
brations can be expressed through the F by means of 



rad 



= — Im 
m 



■ rad 



it) 



N. 



(6) 



where N = Vue is the number of free electrons in the 
MN and V refers to the particle volume. It is necessary 
to underline here that the radiation damping rate is pro- 
portional to the total number of oscillating electrons in 
MN. 

Supposing d(t) = dp exp (— iwi), we obtain for j-th 
component of a radiative linewidth the following expres- 
sion: 



j,rad — „ 



3 mc^ 



Im [1 + [ejj/era - I)] ■ (7) 



Putting here e"^ = we get known expression, e.g., 
from Ref. [H. 

For the sake of simplicity, we will assume that the di- 
electric matrix has no infiuence on the MN and can be 
characterized by 



j(a;) = const 



(8) 



i.e., the dielectric constant of the surrounding medium is 
assumed to be frequency independent. However, it may 
happen in some actual cases that the dielectric medium 
is strongly absorptive at frequencies below the ujpi. If 
that is the case, then the e"j is strongly dependent on 
frequency and contributes to the attenuation of the oscil- 
lations. With accounting for the dielectric matrix prop- 
erties given by Eq. ([5]), Eq. (O can be rewritten as 



j.rad 



6 mc^ 



(9) 



If one will consider the case of the frequency close to 
the frequency of the bulk plasma oscillations of electrons 
in metal Wpi, then the imaginary part of the dielectric 
function tensor for free electron gas can be expressed 
within the Drude-Sommerfeld model aa^^^ 



C7',(U.) 



(10) 



where a'^^ is the principal components of the real part of 



the conductivity tensor. Taking into account the expres- 
sion ((To]), we get 



ri,rad(w) 



3 mc3 ^ 



(11) 



The real part of the conductivity tensor can be ex- 
pressed through the imaginary part of the polarizability 
tensor ajj by means of 

^jj(^) = + - l)l'lm«,,(^). (12) 

In the case of the SPR, one can use in Eq. ([3]) for 
nonspherical MNs, the frequency 



^Soo + [l/L, - l)r 



(13) 
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Here n = Cm and = 1 + einter is the high frequency 
dielectric constant due to interband and core transitions 
of the inner electrons in a MN's material. 

There are different possibilities to calculate a'jj (w) for 
the different frequency regime. Below, we will demon- 
strate how to calculate the a'jj as applied to the nanopar- 
ticle case. 



respectively/24 The semiaxes are connected to the radius 
of sphere R of an equivalent volume through the relation 

The last summand in Eq. ()15p represents the oscilla- 
tion part of the ^ function and the first one refers to its 
smooth part. The function in Eq. p3|) varies with the 
angle because the parameter q becomes dependent on 
the angle 9 for a spheroidal particle, namely 



IV. CALCULATION OF a^^ 

To calculate the tensor a'jj we will use the kinetic equa- 
tions method. Benefit of this method is that it permits 
one to study the effect of the particle shape on the mea- 
sured physical values. Second, it enables us to investigate 
the particles whose sizes are those that the particle sur- 
face start to play an important role. A diffuse boundary 
scattering is assumed to be a good approximation in this 
case. 

The simplest form of a nonspherical shape is a 
spheroid. We restrict ourselves to the nanoparticles with 
a spheroidal shape only. Applying the mentioned method, 
we have found^^ that the components of the conductiv- 
ity tensor for light polarized along (||) or across (_L) the 
rotation axis of a spheroidal MN are 



9uj. 



pi 



Re 



1 



tt/2 



f sin cos^ I 



i sin^ ( 



*(6l) de 



V — Vp 

where v is the electron collision frequency and 9 is the 
angle between rotation axes of the spheroid and direc- 
tion of an electron velocity. Here and below, the upper 
(lower) symbol in the parentheses on the left-hand side of 
Eq. (I14p corresponds to the upper (lower) expression in 
the parentheses on the right-hand side of this equation. 
The subscript v = vp means that the electron velocity 
in the final expressions should be taken on the Fermi 
surface. 

The complex 5* function entering in Eq. (jl4p has the 
form 



4 / 1 

*(<?) = m - ^ ( 1 



where 



$(g) 
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t 1' 



(15) 



(16) 



One can see from Eq. dTG]) that the q is governed by 
the "deformed" electron velocity, which in the case of a 
spheroidal MN takes the form 



v' = vR\ 



sm ( 



cost 
Rii 



v\9), (17) 



where v refers to the electron velocity in a spherical par- 
ticle with radius _R; _R|| and are the spheroid semi- 
axes directed along and across the spheroid rotation axis. 



Vp 



(18) 



A. Conductivity tensor in high-frequency limit 

Let us introduce the frequency of electron oscillations 
between particle walls as 



2R' 



(19) 



Depending on sizes of MN, its shape and temperature, 
the variety of relations between frequencies Ug , v and 
can be achieved. For example, for the Na nanoparticle 
with the radius of i? < 2 A, ~ cjpi. On the other 
hand, with R > 126 A, the electron oscillation frequency 
becomes < where v ~ 4.24- 10^'^ s~^ as is estimated 
, for the Na at SOO'^ K. This leads to different expressions 
for cr(a;), which can be used in Eq. (fTT|) for calculation of 
the plasmon linewidth. 

Below, we will consider the case when the contribution 
of the bulk damping to the radiative plasmon linewidth of 
SPR is neglected. The components of the conductivity 
tensor for a spheroidal MN in the highfrequency (HE) 
limit {lo ^ Vg) and Vs ^ v, can be represented as^^ 



(20) 



where i?j_(= Rx^/^,x = R±/R\\) is a spheroid semi- 
axis directed across to the spheroid rotation axis, and 
riicp) and p{ep) are some smooth functions^^ dependent 
only on the spheroid eccentricity Cp = Vl — a;^ (a prolate 
spheroid), or Cp = y/x^ — 1 (an oblate one). 

One can use the following asymptotic expressions for 
functions ri{x) and p{x) in the cases of both the extremely 
small or the large axial ratio, respectively: 



r]{x) 



p{x) 



tt/8 + 3ttx^ /16, for a prolate spheroid 
x/2 + 1/(42;), for an oblate spheroid 
37r/16 -I- TTx'^ /32, for a prolate spheroid 



X I -1+4 In 2a 
4 ^ 8x ' 



for an oblate spheroid 



(21) 

If one consider the nanowires and nanorods, which can 
be reasonably approximated as prolate spheroids, then 
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one can put r] ~ tt/8 and p ~ 37r/16 with a sufficient 
degree of accuracy. In tlie case of MNs witlr a splrerical 
sliape r/ ~ p ~ 2/3. 



Radiative damping for spheroidal MN in HF limit 



In tlie liiglrfrequency limit, wlren Eq. (j20p can be ap- 
plied, for two components of the linewidth of a spheroidal 
MN embedded in a medium with Cm, we obtain the fol- 
lowing equation 



(i),rad,sp 4^mUJC^RA 



■.L( 



(22) 



where explicit expressions for ri{ep) and p(ep) can be 
found in Ref. IsgI . Accounting for 



(23) 



Eq. ((22)) at the resonance frequency ([T3| takes the form 



r 



pi 



(|),rad,sp 4^3 



(1) 



2n2 



'Kep) 
p(ep) 



(24) 

The depolarization coefficients -^j^ii) for spheroidal MN 

can be found elsewhere, In the case of MNs with a 
very elongated shape {x <^ 1), they can be expressed as 









Hl- 


f) 


- 1 









i±(x) = [l-L||(a;)]/2. 

(25) 



2. Radiative damping for a spherical MN in HF limit 



B. Environment effect 

Because the effect of an electric field on the embedded 
nanoparticles becomes weaker in a dielectric media pro- 
portionally to its refractive index, the environment effect 
plays an important role. The spectral peculiarities of an 
environment effect recently were investigated for the Ag 
and Au nanoparticles, for instance, in Refs. (26ii38i] . 

In order to study the significance of the radiative 
linewidth behavior in more general situations as above 
presented, it is necessary to perform the numerical cal- 
culation in Eq. (fTTj) with the use of a general expression 
for the conductivity tensor given by Eq. (ITil) . Then, the 
radiative linewidth of SPR with an account for Eqs. (He 
and (j23p can be expressed in the form 



(27) 



It allows to account an other important factor, namely 
the dependence of the radiative damping on the particle 
shape. As already was outlined)^'^'^ any change of the 
nanoparticle shape from a sphere, that introducing of 
an anisotropy, results in the splitting of the SPR into 
two modes: a transverse one {T±, perpendicular to the 
spheroid axis of a revolution) and a longitudinal one (Fn , 
parallel to this axis). 

In Fig. 1 we illustrate the behavior of the radiative 
damping components as a function of the medium re- 
fractive index for a fixed particle axes ratio. 

The calculations were conducted for the Au nanopar- 
ticle with the use of Eqs. and dHl), and the fol- 
lowing parameters^: Ue — 5.9 x 10^^ cm"'^, vp = 
1.39 X 10^ cm/s, Wpi = 1.37 x lO^^ s^^ and e^o = 9.84. 



In the case of MN with a spherical shape R\\ = Ri_ = 
R, and one can put the depolarization factor equal to 
L|| = L_L = 1/3 in Eqs. ©, ©, ©-([131), (ED, and 

Eq. (|24p for MN with a spherical shape can be rewrit- 
ten as 

r.ad.p = ^(— j vpR \]^^- (26) 

The increase in linewidth from the radiative damping 
(as we can see from Eqs. ([M)) and ([Ml)) is proportional to 
the MN surface area. This can be understood from the 
fact that the surface scattering is a solely mechanism for 
the change of the electron velocity and acceleration. We 
do not account any other radiative mechanisms in our 
theory as, for instance, the electron-electron collisions. 

The estimations of F for the spherical Au (eoo = 9.84) 
and Na (eoo = 1.14)3- particles with 2R = 400A embed- 
ded in the vacuum (n = 1), in accordance with Eq. ([^5]). 
give: 6.7 and 0.795 meV (or 98 and 1209 fs), respectively. 

If the MN is embedded in the dielectric media with 
Em > 1, then the environment effect ought to be taken 
into account. 
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FIG. 1. (Color online) The radiative linewidth of the surface plas- 
mon resonance components (longitudinal || and transverse ±) vs 
the medium refractive index for a prolate Au nanoparticles with 
the axes ratio iJx/^|| = 0.618. The dashed line corresponds to the 
spherical Au particle with the radius R = 200A. The inset shows 
the same dependence for two spherical Na nanoparticles with the 
radii of 180 and 200A (solid lines). The dashed lines represent 
the r(n) in the case, when the oscillation terms in Eq. II30I I are 
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neglected. 

We observe from Fig. 1 that both radiative damping 
components are decreased with increasing the medium 
refractive index, in other words, the effect is weaker for 
higher dielectric constant of environment. The more of 
MN radius or the less the refractive index is, the higher 
is the calculated radiative damping. 

In general, the resonance plasmon damping in the pro- 
late An nanoparticle was found to be weaker along the 
spheroid revolution axis than the one across this axis. For 
the oblate An nanoparticle, on the contrary, the damping 
along the revolution axis was stronger than that across 
this axis. This result holds regardless of whether the 
photo-excitation is close to the SPR or far from it. 

The resonance peak position is shifted with changing 
the refractive index of the surrounding medium. The 
spectral direction of the shift depends on a number of 
factors, well studied in earlier publications 



1. Radiative plasmon linewidth for a spherical MN in a 
more general case 



The dashed (smooth) curves in the inset in Fig. 1 cor- 
respond to the simple case when both the second and 
the last terms are neglected in the square brackets of 
Eq. (15(1)) . In this simple case, the radiative linewidth of 
the SPR is decreased with enhancing of a refractive index 
n of the environment. This takes place because the ra- 
diative plasmon decay is weakened in the media with the 
higher refractive index owing to decrease of an external 
EM field inside the MN. 

As one can see from the inset in Fig. 1, the oscillations 
of r around its smooth curve make significant correc- 
tions to the smooth picture, especially at large n. The 
oscillations are well pronounced for the Na nanoparticles 
with the small radii and are disappeared for NP with 
a larger radii. This can be connected with the number 
of an electron oscillations between particle walls, which 
is decreased as the particle radius is increased. These 
oscillations is damped markedly with n decreasing and 
practically disappeared at n < 2. The damping enhances 
as the radius of MN becomes larger. 



We consider here the case v' — v, when the 5* function 
ceases to depend on the angle 6. Then the integration 
over 6 in Eq. gives 1/3 and we derive for a' 



(28) 



with q = 2R{v — iiS)lv taken at w = vp. 

Let us choose the case v ^ Vs^ for illustration. After 
simple algebraic calculations with the use of Eqs. ((T5]) 
and II16D. we obtain in this case: 



2vs UJ Iv'- I UJ 
1 sm 1 1- I 1 — cos — 



Substituting Eqs. ([29|) and into Eq. ((TTl), we get 



(29) 



sph 



9< V C 



2 2 
1- -sinC+— (1-cosO 



with 



e = c(",i?) = 



2i?a;pi 



(30) 
(31) 



One can see from Eq. (j30| that the radiative damping 
in a spherical MN oscillates with altering of both the 
radius of MN and/or the refractive index of an embedding 
medium. Eq. (pO| at ^ ^ 1 transforms into Eq. (|26)) . 

The inset at the Fig. 1 shows the behavior of the SPR 
linewidth for spherical Na nanoparticles with different 
radii. For calculation we use Eq. (1501) and the following 
parameters: ~ 2.65xl02'^ cm^'^, vp = 1.07x10^ cm/s, 
Wpi = 9.18 X 10^^ s-\ and too = 1-14. 



2. Ejfect of the particle size 

Figure 2 depicts the behavior of the radiative linewidth 
of SPR as the function of MN radius for some fixed 
medium refractive indices. The calculations were fulfilled 
for the Au nanoparticle with the use of Eq. ([50]) and the 
same numerical parameters as given above. 

As we can see in Fig. 2 and from Eq. (PO)) . the F is 
increased quadratically with R (dashed lines). We sup- 
pose, this is due to the fact that the radiative damping in 
the MN with a spherical shape is proportional to the area 
of a sphere. The growth in F occurs slower in the media 
with a higher refractive index. The redaction of radiat- 
ing damping in the media with higher n implies a reduced 
dephasing of the plasmon mode. Because the effect of an 
external electric field on the embedded nanoparticle be- 
comes weaker in a dielectric media with higher refractive 
index, the resonance radiative damping tends to decrease 
within the nanoparticle. 

The F oscillates around a smooth curves with an in- 
crease of the particle radius. The period of these oscilla- 
tions is enhanced for the MNs embedded in the medium 
with a higher refractive index. The magnitude of these 
oscillations is the greater the higher refractive index is. 
The oscillations of F with R follow the quadratic depen- 
dence as well. The inset at the Fig. 2 shows the same 
dependence for Na nanoparticle embedded in the media 
with n — 1 and n = 9. 

On the oscillations of the SP lifetime as a function of 
nanoparticle size it was pointed out for the first time 
in Ref. ,21, where the semiclassical theory was used to 
evaluate the SPR in MNs. 
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FIG. 2. (Color online) The radiative linewidth of SPR vs radius of 
spherical Au nanoparticle (in units of the Bohr radius ag = 0.53) 
embedded in the water n ~ 1.33 or in the medium with refractive 
index n = 9. The dashed lines represent the r(i?) in the case, when 
the oscillation terms in Eq. I I3UII are neglected. The inset shows the 
same dependence for Na nanoparticle. 

The behavior of T at small particle radii R/as < 40 
in our approach practically does not depend on both the 
particle radius and the medium refractive index and can 
be described more precisely in the frame of the quantum 
theory. 

So far we have considered only radiative processes. Be- 
low, we dwell shortly on the nonradiative processes when 
the electron scattering in the MN dissipates oscillation 
energy into heat. 



So, the decay time of the plasmon resonance is the opti- 
cal conductivity of the MN at the light frequency multi- 
plied by a geometrical factor. To study the nonradiative 
linewidth for MNs with a given Lj , it is enough to calcu- 
late the real part of the conductivity tensor as a function 
of frequency. 

As one can see from Eq. ([32]) . the radiative damping 
rate start to dominate as the light scattering cross sec- 
tion by MN exceeds the light absorption one in it. The 
calculations for spherical MN in vacuum gives that radius 
of MN for which both cross sections become comparable 
each with other. In the high-frequency limit, we found 
that 

R„,^(i)"\c.,r"^. (35) 

For instance, R ~ lOOA for Au particle at the plasmon 
frequency uj = Wpi / \/3 . 

In the case of low frequencies, we have 

Rlfc^—.^^. (36) 

The latter formula gives R ~ 7860A for Au particle em- 
bedded in a medium with = 1. 

In order to take into account the radiative damping to- 
gether with collisions of free carriers with the MN surface, 
the effective damping rate TeS = rnonrad + Fjad must be 
introduced. For understanding of the decay mechanism 
of the electron plasma oscillations the knowledge of the 
decay time is of central importance. 



C. Linewidth for a nonradiative processes 

Another way for calculation of Fj iad supposes the 
knowledge of the value of a nonradiative damping 

-^j,nonrad- That is 



- j,rad 



CTabs 



- j,nonrad7 



(32) 



where Csca and Cabs are the light scattering and absorp- 
tion cross sections, respectively. 

The nonradiative damping rate of a SP can be ex- 
pressed as^ 



rj,nonrad(^) 



At: Li 



<,(w), 



(33) 



with the a'jj given in the Section IV. Eq. ()33|) defines 
the linewidth or, correspondingly, the decay time of the 
plasmon resonance due to electron scattering both from 
the bulk and from surfaces of the particle. 

In the case of medium with e'„ — >• 1, Eq. (|33l) is reduced 
merely to 



rj,nonrad(w) = AttLj (t'^Auj) 



(34) 



V. SUMMARY AND CONCLUSIONS 

We used a local field approach and a kinetic equation 
method to study the plasmon resonance linewidth for 
metal nonspherical nanoparticles embedded in any dielec- 
tric media. It enables to calculate the radiative linewidth 
for MNs with different geometry with an account for the 
light scattering from the particle surfaces. 

The general formula is proposed for a damping rate or 
a decay time due to electron scattering from the bulk and 
particle surfaces. By means of this formula one can eval- 
uate the linewidth directly through the tensor of optical 
conductivity of the MN. 

With changing the MNs shape from spherical to the 
spheroidal one, the single radiative plasmon resonance 
splits into two components: the longitudinal and a trans- 
verse one to the spheroid rotation axis. We found that 
both the radiative damping components for Au nanopar- 
ticle are decreased with increasing the medium refractive 
index. The resonance plasmon damping in the prolate Au 
nanoparticle was found to be weaker along the spheroid 
revolution axis than that across this axis. 

For Na nanoparticle, we detect the oscillations of F 
with the refractive index increasing. The amplitude of 
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these oscillations enhances for media with higher dielec- 
tric constant. The oscillations arc well pronounced for Na 
nanoparticles with the small radii and are disappeared for 
Na nanoparticles with a larger radii. 

We clearly show for spherical MNs that the radiative 
linewidth of SPR enhances quadratically with the parti- 
cle radius increasing. The growth in T occurs slower in 
the media with a higher refractive index. The T oscillates 
as well with changing in the particle radius. The mag- 
nitude of these oscillations is enhanced markedly for the 
MNs embedded in the medium with a higher refractive 



index. 

The effects of both the particle radius and the environ- 
ment on the radiative plasmon resonance linewidth are 
illustrated by the example of Au and Na nanoparticles. 

The contribution of the nonradiative plasmon decay is 
discussed as well. 

Our theoretical results should be important for the 

analysis of the experimental data on the optical and 
transport properties of MNs embedded in various dielec- 
tric media. 
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